Mechanism for Odd Parity Superconductivity in Iron-Based Superconductors 
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Under the assumption that superconducting pairing is driven by local d-p hybridization, we show 
that the superconducting state in iron-based superconductors is classified as an odd parity s-wave 
spin-singlet pairing state in a single trilayer FeAs/Se, the building block of the materials. In a low 
energy effective model with only d-orbitals in an iron square bipartite lattice, the superconducting 
order parameter in this state is a combination of a s-wave normal pairing between two sublattices 
and a s-wave ^-pairing within the sublattices. Parity conservation was violated in proposed super- 
conducting states in the past. The results demonstrate iron-based superconductors being a new 
quantum state of matter and suggest that a measurement of odd parity can establish fundamental 
principles related to high temperature superconducting mechanism. 



In a strongly correlated electron system, major physics 
is determined locally in real space. Important properties, 
such as pairing symmetry in a superconducting state, are 
expected to be robust against small variation of Fermi 
surfaces in reciprocal space. Although superconducting 
mechanism related to high temperature superconductors 
(high T c ) is still unsettled, the robust d-wave pairing 
symmetry in cupratespQ can be understood under this 
principle. 

Is this principle still held for iron-based 
superconductors [2H5]? Namely, do all iron-based 
superconductors possess one universal pairing state? 
Unlike cuprates, the answer to this question is highly 
controversial because different theoretical approaches 
have provided different answers and no universal state 
has been identified [6 . Nevertheless, as local electronic 
structures in all families of iron-based superconductors 
are almost identical and phase diagrams are smooth 
against doping [5-7 , it is hard to argue that the materials 
can approach many different superconducting ground 
states. 

In a recent paper [8 , one of us provided a complete 
symmetry classification for pairing symmetries in iron- 
based superconductors and showed that spin-singlet pair- 
ing with odd parity can naturally take place because of 
the intrinsic 2-Fe unit cell, which was ignored in pre- 
vious theoretical studies based on effective models with 
1-Fe unit cell. It was further argued that the pairing 
between electron pockets is most likely Ai u (D 2 d) s-wave 
or £?2u {Cav) d-wave 77-pairing. The pairing state has 
both s-wave characters, such as no symmetry-protected 
nodes or nodal lines on superconducting gap functions, 
and d-wave characters, such as a sign change between 
top and bottom As/Se layers in real space. However, no 
microscopic mechanism has been proposed for odd parity 
pairing. 

In this Letter, we provide a microscopic understand- 
ing to show that the superconducting state in iron-based 
superconductors is classified as an odd parity s-wave spin- 
singlet pairing state (OPS) in a single trilayer FeAs/Se, 



the building block of the materials. This conclusion is 
only based on the consensus that As(Se) plays a critical 
role in driving superconductivity so that superconducting 
pairing is determined locally by d-p hybridization. Under 
this consensus, in a low energy effective model with only 
d-orbitals in an iron square lattice, we show that a nor- 
mal pairing order defined between two iron sublattices 
is parity odd. Namely, it carries a sign change between 
top and bottom As/Se layers. This is because the effec- 
tive nearest neighbor (NN) hopping of t^g d-orbitals in- 
duced by d-p hybridization is generated through the anti- 
bonding state of p-orbitals formed between top and bot- 
tom As/Se layers. As parity is a good quantum number, 
a superconducting order parameter between two next NN 
(NNN) sites must be an ^-pairing. Thus, in the effective 
d-orbital models, a superconducting state that does not 
violate parity conservation must include both normal and 
r] pairing. An OPS, classified as a Ai u {D 2 d) s-wave or 
B2u (C± v ) d-wave in a full lattice symmetry [8], is a combi- 
nation of a s-wave normal pairing between two sublattices 
and a s-wave 77-pairing within the sublattices. We show 
that the meanfield Hamiltonian of this state explains the 
dual s-wave and d-wave type characters observed exper- 
imentally and unifies the description of iron-pnictides[7 
and iron-chalcognides[9HT2] . The results conclude that 
parity conservation was violated in proposed supercon- 
ducting states [6] in the past. The confirmation of the 
odd parity state will have a tremendous impact on un- 
derstanding high-T c mechanism. 

Before we start our analysis, we repeat the definition 
of normal and rj pairing in [8]. Let k be momentum with 
respect to the 1-Fe unit cell in an iron square lattice. The 
normal pairing refers to (—k,k) pairing and the 77-pairing 
refers to (—k,k + Q) pairing where the momentum vector 
Q — (7r, 7r) is a reciprocal lattice vector in a 2-Fe unit cell. 

Effective Hamiltonian We consider a general Hamilto- 
nian in a single trilayer Fe — As(Se) structure coordi- 
nated by Fe and As(Se) atoms, 

H = Hdd + Hdp + H pp + Hj (1) 
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FIG. 2: (a) The calculated band structure of FeSe with 
the weight of p orbitals of Se. (b) The decomposed charge 
density of the band at Y marked by letter B (anti-bonding 
states). (c)The decomposed charge density of the band at M 
marked by letter C (bonding states). 



FIG. 1: (Color online) The nearest neighbor hopping param- 
eters for intra-dry-orbital are shown in (a) and (b). tdd is the 
amplitude of the direct hopping of d xy orbital while tdpd is 
the amplitude of the indirect hopping through p z orbital of 
As/Se atom. The reason of sign change for td P d between (a) 
and (b) is that the p z orbitals in top layer and bottom layer 
form occupied bonding states in (a) and empty anti-bonding 
states in (b) . The difference is illustrated by filled and empty 
p z orbitals in (a) and (b). (c) and (d) shows the local p — d 
s-wave pairing pattern and the gauge transformation between 
them. 



where Hdd, Hd p and H pp describe the direct hopping be- 
tween two d-orbitals, the d — p hybridization between 
Fe and As(Se) and the direct hopping between two p- 
orbitals respectively. Hi describes any standard interac- 
tions. Here we do not need to specify the detailed param- 
eters. This Hamiltonian has a full symmetry defined by 
a non-symmorphic group which can be specified equiva- 
lently as G = {E, I) ® C± v or G = (£", I) ® D 2 d as shown 
in [8] , where / is the space inversion operation defined at 
the center of a NN Fe link, D 2 d is the point group at iron 
sites and C^ v is the point group at the center of an iron 
square. 

An effective Hamiltonian is obtained by integrating out 
p-orbitals, which can be written as 



eff 



dd,eff 



Hj 



,eff- 



(2) 



The effective band structure can be written as Hdd,eff = 
Hdd + Hdpd, where Hd P d is the effective hopping induced 
through d-p hybridization. Hdd,eff has been obtained by 
many groups [ 13HT7] . The major effective hopping terms 
in Hdpd can be divided into two parts Hdpd,NN, which 
describes NN hopping and Hd P d,NNN , which describes 
NNN hoppings in the iron square lattice. If one carefully 
checks the effective hopping parameters for t 2g orbitals in 
Hdpd,NN, one finds that they have opposite sign to what 
we normally expect in a natural gauge setting as shown in 
figjTJa,b), where d xy orbital is illustrated as an example. 



We see that the hopping parameter tdd must be negative. 
However the effective hopping parameter, td p d, is posi- 
tive and even larger than \tdd\ in [T5HT7] . In a tetragonal 
lattice, tdpd can only be generated through d xy — p z hy- 
bridization. A positive value of tdpd suggests that virtual 
hopping which generates tdpd must go through an unoc- 
cupied p z state. As shown in figjTJa,b), a d xy equally 
couples to p z orbitals of top and bottom As atoms. A 
high energy p z state must be an anti-bonding p z state 
between NN As atoms. This analysis is held for all t 2g 
orbitals which play the dominating role in low energy 
physics. It is also easy to check that the effective NNN 
hoppings between t 2g orbitals are dominated through an 
occupied p states, which is primarily a bonding state of p 
-orbitals. Therefore, the NN effective hoppings are gener- 
ated through d—p a hybridization, where p a represents an 
anti-bonding p-orbital states and the NNN effective hop- 
pings are generated through d — p^ hybridization where 
Pb is the bonding p-state. 

The above microscopic understanding is not surpris- 
ing. In fact, it is known in LDA calculations [T4| fT8| [T9] 
that p-orbitals in As/Se are not fully occupied and there 
are significant overlappings between p-orbitals on bot- 
tom and top As/Se layers. Moreover, since Hd P d,NN and 
Hdpd,NNN primarily affect hole pockets around Y and 
electron pockets at M separately, we can check the dis- 
tribution of anti-bonding p states and bonding p states 
in band structure to further confirm the analysis. In 
figj2ja), we plot the band structure of FeSe and the dis- 
tribution of p orbitals. The p z orbitals of Se are mainly 
at +1.5 eV at T and -3 eV at M. By analyzing the bands 
at T and M as shown in Fig(2^b) and (c), we confirm that 
the p z orbitals of Se at Y and M belong to anti-bonding 
and bonding states separately. 

Hidden Z 2 symmetry characters in effective Hamilto- 
nian: Knowing the above hidden microscopic origins in a 
derivation of an effective Hamiltonian allows us to under- 
stand the symmetry characters of the effective Hamilto- 



3 



nian in the original lattice symmetry. As shown in|5], in 
the original lattice symmetry, G, we can introduce a Z 2 
classification specified by a^, where &h is the reflection 
along z-axis. 

The original Hamiltonian is invariant under (<j^,T), 
where T is an in-plane translation by one Fe-Fe lattice. 
The d — p a hybridization is odd under &h while the d — p^ 
hybridization is even under &h- Thus, the NN hopping 
Hd P d,NN and NNN hopping Hdpd,NNN should be classi- 
fied as odd and even under &h respectively. Namely, 

&hHdpd,NNO'h = ~ 1 

&hHdpd,NNN&h = 1 (3) 

The above hidden symmetry property is against the main 
assumption taken in many weak coupling approaches, 
which assume that the essential physics is driven by the 
interplay between hole pockets at T and electron pockets 
at M [6 j. As indicated in figj2^a), the interplay between 
the hole and electron pockets must be minimal because 
of their distinct microscopic origins. 

Gauge transformation and parity of pairing order pa- 
rameter: The symmetry difference in eq(3]has fundamen- 
tal impact on how to consider the parity of a supercon- 
ducting state if superconducting pairing is driven by local 
d-p hybridization. 

It has been shown that in a system where short range 
pairings in real space dominate, superconducting order 
parameters are momentum dependent and a gauge prin- 
ciple must be satisfied because the phases of supercon- 
ducting order parameters can be exchanged with those 
of the local hopping parameters [20} |2T] by gauge trans- 
formations. As an example, a d-wave superconduct- 
ing state in cuprates can be mapped to a s-wave su- 
perconducting state by a gauge mapping which changes 
the hopping terms from s-type symmetry to d-type 
symmetry [22 . Therefore, only the combined symmetry 
of hopping terms and pairing orders associated to them 
is a gauge-independent symmetry character to classify 
states. This gauge principle does not exist in a conven- 
tional BCS-type superconductor in which the informa- 
tion of pairing in real space is irrelevant. 

Now we apply the gauge principle and let Ann and 
Annn be superconducting order operators associated 
with Hdpd,NN and Hdpd,NNN respectively. In a super- 
conducting state, we must have 

[A AT at] [Hd p d,NN] = [ANNN][Hdpd,NNN] (4) 

where [A] indicate the symmetry of A. Following eqjij 
we have [A at at] = —[A at at at] under &h- Therefore, based 
on the classification of pairing symmetries in [8], we im- 
mediately conclude that the parity is odd and the su- 
perconducting order < A at at at > must be an 77-pairing if 
< Ann > is a normal pairing. 

The above analysis can be easily illustrated in real 
space. As shown in fig|3j if superconducting pairing is 




FIG. 3: (Color online) The NN and NNN p-d local pairing 
patterns with odd parity are shown in (a) and (b) in the 
natural gauge. Note that p orbit als of As/Se in (a) form the 
anti-bonding states while that in (b) form the bonding states. 
We distinguish the two states with different filled green balls 
between (a) and (b). The p — d pairings can be projected into 
effective d — d pairings shown in the bottom row. 

driven by local d-p hybridization, the superconducting 
order is a pairing between d and p orbitals Ad p =< 
d + p + >. A uniform < d + p+ > is parity odd. The NN 
pairing, < Ann > in the effective model, must originate 
from < d + p+ > and thus is also parity odd. The gauge 
principle can be understood as shown in fig(TJc,d). If we 
can take a new gauge for Fermion operators of p-orbitals, 
p —> — p, in one of the two As(Se) layers, the anti-bonding 
operator p a maps to the bonding operator p^. This gauge 
mapping exactly transfers the parity between hopping 
terms and superconducting order parameters. 

Meanfield Hamiltonian for odd-parity s-wave state: 
The above analysis can be generalized to all effective hop- 
pings. The basic idea is to divide the iron square lattice 
into two sublattices. In an odd parity s-wave state, the 
pairing between two sublattices must be normal pairing 
while the pairing within sublattices must be ^-pairing. 
Therefore, a meanfield Hamiltonian to describe the odd 
parity s-wave state in 1-Fe unit cell can be generally writ- 
ten as 

H m f = H d d,eff + ^2 (da(3,nA a p jn (k) + 5 a p iV ^a/3,r] (k) + ft.C.)(5) 

a,P,k 

where A a ^ n = d a ^(k)d^(-k) - d ai (k)d^(-k) and 
Aaf3, v = d at (k)df3i(-k + Q) - d a i(k)dp^(-k + Q). In 
general, the normal and r] pairing order parameters sat- 
isfy 

<W,n(>) = -<W,n(& + Q) (6) 

Sapvik) = S a p, v (k + Q). (7) 
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FIG. 4: (Color online) The Fermi surfaces of a five-orbital 
model in [15] are shown in (a) and (b). The forms of hopping 
terms and hopping parameters can be found in [15]. Here, we 
only add a chemical potentials to tune the fermi level. We 
set \i — 0.1 and —0.04 in (a) and (b). The high-symmetry 
points are shown in (a) and the pairing channels connect- 
ing the points on the fermi surface are denoted by the black 
lines with arrows. The quasi-particle spectrum of the super- 
conductive states for (a) and (b) are shown in (c) and (d). 
We can find the (c) is full gaped and (d) has nodes at the 
electron pockets. The superconductive order parameters are 
chosen: Aft x = Aft y = 0.05; Aft 4 = 0.05; Aft^ = 0.05; 
Aft" = 0.05;^ = -0.1; 



These equations capture the sign change of supercon- 
ducting order parameters in momentum space. The sign 
change here is required by symmetry. As the inter- 
orbital pairing can be ignored for s-wave pairing and 
the pairing is dominated by NN and NNN pairings, the 
important parameters are 5 aa , n oc cosk x + cosk y and 
5aa,r] oc cosk x cosk y . Thus, the superconducting gaps on 
hole pockets are mainly determined by S aajU and those on 
electron pockets are mainly determined by 5 aa ^. While 
detailed studies will be carried out in the future, the su- 
perconducting gaps obtained from eq|5] can explain ex- 
perimental results observed in both iron-pnictides and 
iron-chalcogenides[23 . There is no symmetry protected 
node in this superconducting state. However, accidental 
nodes can easily take place. In fig(4j we plot numerical 
results for two cases. Parameters are specified in the cap- 
tion of the figure. The superconducting gap in the first 
case is a full gap while it has gapless nodes on electron 
pockets in the second case. 

Impact to High-T c Mechanism In ref.jS], the odd par- 
ity s-wave state was conjectured based on intriguing ex- 
perimental facts. With microscopic mechanism proposed 
here, we can address important impact on high T c mech- 
anism for iron-based superconductors and other high T c 
superconductors if the state is confirmed. 



First, the microscopic mechanism revealed here fun- 
damentally differs from those proposed in weak coupling 
approaches which only emphasize Fermi surfaces. Fermi 
surfaces are only determined by energy dispersion. It 
provides no information about underlining microscopic 
processes which are local and bound with high energy 
physics. In correlated electron systems, these processes 
essentially determine many important properties. 

Second, our study provides fundamental reasons why 
we failed to recognize the odd parity symmetry in the 
past. The symmetry principle and gauge principle were 
not properly handled. In the past, the effective Hamil- 
tonian was viewed in the symmetry g roup C4.1) at iron 
sites rather than the original lattice symmetry G. We 
can see that if &h could be set to one, G is reduced to 
C^y. However, due to the anti-bonding p orbital states, 
the effective Hamiltonian does not represent correct sym- 
metry of the original lattice in a natural gauge setting. 
The correct physics can only be understood after the hid- 
den gauge is revealed. For an order parameter which is 
momentum dependent, this gauge information is critical. 
The gauge principle becomes very important for us to 
search new physics in other complex electron systems. 

Finally, if the odd parity state is confirmed, the funda- 
mental objects in superconducting states of high T c ma- 
terials must be the tightly binding Cooper pairs between 
d and p orbitals. In this view, the odd parity s-wave 
state closely resembles the d-wave state in a Cu-0 plane 
of cuprates. We expect an identical mechanism to select 
sign changed superconducting orders in both materials. 

In summary, we provide a microscopic mechanism to 
support the formation of an odd parity s-wave supercon- 
ducting state in iron-based superconductors. We demon- 
strate that in an effective model based on d-orbitals, both 
normal pairing and n pairing must be included if the su- 
perconducting state conserves parity. Superconducting 
states studied in the past violate parity conservation. 
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